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I. INTRODUCTION 

Interest in the field of charmonium spectroscopy has 
been renewed with the recent discovery of numerous 
charmonium and charmonium- like states fH-flll]. How- 
ever very little is known about the first radially ex- 
cited x c j(2P) states which are expected to exist in the 
mass region from 3.9 to 4.0GcV/c 2 , just above the DD 
threshold [l2j . The Belle collaboration has observed the 
Z(3930) state in 77 production of the DD system [l]|, 
and this is considered a strong candidate for the Xc2(2P) 
state; indeed it is so- labelled in Ref. [T3|. The Belle 
analysis obtained a mass of to = (3929 ± 5 ± 2) McV/c 2 
and a total width of T = (29 ± 10 ± 2) MeV, with 
quantum numbers J PC = 2 ++ preferred. The par- 
tial width T 77 x S(Z(3930) — > DD) was determined as 
(0.18 ± 0.05 ± 0.03) keV, where T 77 is the radiative width 
of the Z(3930) state, under the assumption that J = 2. 
The observation of this state has not been confirmed so 
far 0. 

In this paper the process 77 — > DD, illustrated by 
the Feynman diagram shown in Fig. [1] is studied in a 
search for the Z(3930) state. In Fig. [TJ the initial state 
positron, e + (electron, e _ ), emits the virtual photon 7]" 
(72), yielding the final state positron, e + ' (electron, e"'); 
the momentum transfer to 7^ (72) is q\ (92)- The virtual 
photons interact to produce the DD final state. When 
the e +/ and the e~' are emitted along the beam directions 
the values of q\ and q\ are predominantly close to zero, 
and the two photons can be considered to be quasi-real. 
Since in this case neither the e +/ nor the e~' are detected, 
the analysis is termed untagged. 

The data sample used in this analysis corresponds 
to an integrated luminosity of 384 fb -1 recorded at the 
T(4S*) resonance (10.58 GeV) and at a center of mass 
(cm.) energy of 10.54 GeV by the BABAR detector at the 
PEP-II asymmetric-energy e + e~ collider. 

The BABAR detector is described briefly in Sec. [JJJ and 
the principal criteria used in the selection of candidate 
two-photon-interaction events are discussed in Sec. IIIII 
The reconstruction of DD pair events is presented in 
Sec. IIV1 and the relevant Monte Carlo simulations are 
detailed in Sec. [V] The purity and reconstruction effi- 
ciency of the 77 — > DD event sample are considered in 
Sees. |VT] and IVII1 respectively, and the DD signal yield 
and invariant mass resolution are presented in Sec. IVIIII 
the mass and total width for the Z(3930) state are ob- 
tained from a fit to the DD invariant mass distribution. 
The angular distribution in the DD rest frame for the 
Z(3930) mass region is studied in Sec. IIX| and the impli- 
cations for the spin of the Z(3930) state are discussed. In 
Sec.|X]thc partial radiative width of the Z(3930) state is 
extracted. Sources of systematic uncertainty are detailed 
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FIG. 1: Two-photon production of the DD system. 

in Sec. lXIl and the results of the analysis are summarized 
in Sec. EH 

II. THE BABAR DETECTOR 

The BABAR detector is described in detail else- 
where [l5|. Charged particles arc detected, and their 
momenta measured, with a combination of five layers 
of double-sided silicon microstrip detectors (SVT) and 
a 40-laycr cylindrical drift chamber (DCH), both coaxial 
with the cryostat of a superconducting solenoidal mag- 
net which produces a magnetic field of 1.5 T. Charged 
particle identification is achieved by measurements of the 
energy loss dE/dx in the tracking devices and by means 
of an internally reflecting, ring-imaging Chcrcnkov de- 
tector (DIRC). Photons and electrons are detected and 
their energies measured with a CsI(Tl) electromagnetic 
calorimeter (EMC), covering 90% of the 4ir solid angle 
in the T(45) rest frame. The instrumented flux return 
of the magnetic field is used to identify muons and . 

III. SELECTION OF 
TWO-PHOTON-INTERACTION EVENTS 

The selection of two-photon-interaction events for an 
untagged analysis is based on established procedures (see 
for instance Refs. 0, [l3)- Due to the small scattering 
angles involved, most of the incoming beam energy is 
carried away by the e +/ and e~' (see Fig. [TJ. This results 
in a large value of the missing mass squared 

TO miss = (Pe+ + Pe- ~ PD ~ PdY ( l ) 

where p e ± are the four-momenta of the beam electron 
and positron and po, Pf) are the four-momenta of the 
final state D and D mesons, respectively. In addition, 
for these events, the resultant transverse momentum of 
the DD system p t (DD) is limited to small values. 

In order to establish selection criteria for 77 — > DD 
events, the reaction 

e+e~ K + K~tt + tt-X (2) 

is studied first using a data sample corresponding to an 
integrated luminosity of 235 fb -1 . The system X con- 
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tains no additional charged particles. This reaction has 
been chosen because it has the same particle configura- 
tion as one of the final states we consider in this analysis. 
The charged kaons and pions are identified as described 
in detail in Sec. lIVI Neutral pions are reconstructed from 
pairs of photons with deposited energy in the EMC larger 
than 100 MeV. It is required that no ir° meson candidate 
be found in a selected event. 

Two-photon production of the K~ K + tt + it~ system 
should yield large values of m x , the missing mass 
squared, 

m X = (Pe+ +Pe- -PK+ ~Pk- -Pit+ -Ptt-Y- (3) 

In addition, production of the K~ K + tt + it~ system via 
Initial State Radiation (ISR) should yield the small val- 
ues of m x associated with the ISR photon, for which 
detection is not required. The observed distribution of 
the K~ K + tt + it~ invariant mass, m(K~ K + tt + t:^ ), re- 
sulting from the reaction of Eq. @ is shown in Fig.^a). 

There are clear signals corresponding to the produc- 
tion of rj c (lS), x c o(lP), and Xc2(l-P), and, since these 
states all have positive C-parity, it is natural to asso- 
ciate them with two-photon production. Similarly, the 
large J/tp signal observed would be expected to result 
from ISR-production, because of the negative C-parity 
of the J/t/j . For the parameters of these states, see Ta- 
ble i 

The distribution of m x for 2.8 < m(K" K+n+n-) < 
3.8 GeV/ c 2 is shown in Fig.^b) . The large peak near zero 
is interpreted as being due mainly to ISR production of 
the K~K + 7r + 7r~ system, while two-photon-production 
events would be expected to occur at larger values of 
m x . This is shown explicitly by the distributions of 
Figs. ^c) and ^d), which correspond to the require- 
ments rn\ < 10 (GeV/c 2 ) 2 and m\ > 10 (GcV/c 2 ) 2 , 
respectively. 

In Fig. [3fc) there is a large J/tjj signal, and a much 
smaller tl'(2S) signal can also be seen. For e + e~ colli- 
sions at a cm. energy 10.58 GeV, the ISR production 
cross section for J/t/j is about three times larger than 
for Tp(2S); also B(J/tp K~ K + tt + h~) is approximately 
nine times larger than the corresponding iJj(2S) branch- 
ing fraction value [13] . 

It follows that the observed J/ip signal would be ex- 
pected to be w 27 times larger than that for ip(2S). The 
signals in Fig.^c) seem to be consistent with this expec- 
tation, and they are also in agreement with the detailed 
analysis of ISR production of the K~ K + ir + ir~ system in 
Ref. [l8[. There is a Xc2(l-P) signal in Fig. |^c) which is 
comparable in size to the ip(2S) signal. The branching 
fraction for tp(2S) K~K+tt+tt- is m 7.5 x 10" 4 QJ], 
while the product B(ip{2S) -> 7Xc2jQ--P)) x B(xtiz{lP) -> 
K-K + TT + n-) is w 7.8 x 10~ 4 0, so that the pres- 
ence of such a Xc2(l-P) signal is consistent with the ex- 
pected transition rates. For the Xci(l-P), B(tp(2S) — > 
lXd{lP)) x B(xci{lP) -► K-K+7T+7T-) « 4.0 x 10" 4 , 
and so a Xci(l-P) signal of approximately half the size of 
the x C 2(lP) signal would be expected in Fig. [5Jc); again 



the data seem to be in reasonable agreement with this 
expectation. 

Finally, for the Xo o(lP), B{^(2S) -> 7Xco(lP)) x 
B(Xco(lP) -> K-K+ir+Ti-) w 16.8 x 10~ 4 , and the cor- 
responding signal in Fig. (2jc) would be expected to be 
about twice the size of the ip(2S) signal. The Xco(l-P) 
signal seems to be larger than that of the ip(2S), but 
not by a factor of two; this may be because the larger 
energy photon from the '4>{2S) — > jxcd(1P) transition, 
when combined with the ISR photon, can yield a value 
of m x which is larger than 10 (GeV/c 2 ) 2 . In summary, 
the signals observed in Fig. [5Jc) appear consistent with 
those expected for an ISR production mechanism, espe- 
cially since there is no indication of any remnant of the 
large r/ c (lS) of Fig. HJa). Furthermore, the XcJ signals 
in Fig. [U[c) arc removed by requiring that the trans- 
verse momentum of the K~ K + ir + ir~ system be less than 
50MeV/c (see discussion of Fig. [UJd) below), which indi- 
cates clearly that they do not result from two-photon 
production. 

In Fig. [U^d), the r] c (lS) signal of Fig. appears to 
have survived the m 2 x > 10 (GcV/c 2 ) 2 requirement in its 
entirety, and the Xco(l-P) an d Xe2(l-P) signals have been 
reduced slightly, as discussed in the previous paragraph; 
in both Fig. (5Ja) and HJd) there is some indication of 
a small signal in the region of the r] c (2S) mass. A J/tp 
signal of about one third of that in Fig. Ufa) is present 
also in Fig. [2jd). This is interpreted as being primar- 
ily due to a) the emission of more than one initial state 
photon, with the consequence that values of m x greater 
than 10 (GeV/c 2 ) 2 are obtained, b) the ISR production 
of the i){2S) with subsequent decay to J/ip + neutrals, 
and c) two-photon-production of the Xc2(l-P) followed 
by Xc2 (If) — > jJ/ip, which has a 20% branching frac- 
tion Q. 

It follows from the above that the requirement m x > 
10 (GeV/c 2 ) 2 significantly reduces ISR contributions to 
the K~ K + ir + ir~ final state while leaving signals associ- 
ated with two-photon-production essentially unaffected. 
For this reason, the requirement that ?n 2 niss of Eq. ([1]) be 
greater than 10 (GeV/c 2 ) 2 is chosen as a principal selec- 
tion criterion for the isolation of events corresponding to 
77 -> DD. 

As mentioned above, it is expected that for an un- 
tagged analysis of 77 — > DD, the transverse momen- 
tum pt(DD) should be small. In order to quantify this 
statement, the data of Fig. [2jd) were divided into in- 
tervals of 50 MeV/c in the transverse momentum of the 
K~ K + ir + ir~ system with respect to the e + e~ collision 
axis, which is considered also to be the collision axis for 
two-photon-production events. For each interval a fit 
was made to the m[K ~K + it + t{~) mass distribution in 
the mass region 2.7 < m{K~ K+tt+tt-) < 3.3 GeV/ c 2 . 
The function used consists of a second-order polynomial 
to describe the background, a Gaussian function for the 
J/tjj signal and a Breit-Wigner for the rj c (lS) signal con- 
volved with a Gaussian to account for the resolution. The 
Pt-dependence of the resulting ?7 C (1S) yield is shown is 
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Fig. [3(a), and that of the J ftp yield is shown in Fig. [3(b). 
The shapes of the distributions are quite similar for 
p t > lOOMeV/c, but the interval from 50 - lOOMeV/c 
contains w 180 more rj c (lS) signal events, and that for 
— 50MeV/c exhibits an excess of ss 800 signal events. 
This behaviour is expected for two-photon-production of 
the T] C (1S). Thus, the requirement p t (DD) < 50MeV/c 
is imposed as the second principal selection criterion for 
the extraction of 77 — > DD events. 

Since the two-photon reactions 77 — > K~ A" + 7r + 7r~ 
and 77 — > DD arc quasi-exclusive in the sense that only 
the final state e + and e~ are undetected it is required 
in both instances that the total energy deposits -Eemc in 
the EMC which are unmatched to any charged-particle 
track be less than 400 MeV. The net effect is a small 
reduction in the smooth background. The histogram of 
Fig. [3jc) corresponds to the K~ K + ir + ir~ candidates of 
Fig. 15(d) after requiring p t (K -R+tt+tt-) < 50MeV/c 
and that the EMC energy sum be less than 400 MeV. 
The p t criterion reduces the r] c (lS) signal by a factor 
ss 2, while the J/tp signal is reduced by a factor w 5, as 
is the continuum background at 2.7GcV/c 2 . More sig- 
nificantly, the continuum background at 3.7GcV/c 2 , just 
below the DD threshold, is reduced by a factor « 10. 

It follows that the net effect of the three principal se- 
lection criteria described above (missing mass rn^ ias > 
10 (GeV/c 2 ) 2 , resultant transverse momentum p t (DD) < 
50MeV/c and total energy deposit in the calorimeter 
Eemc < 400 MeV) is to significantly enhance the number 
of two-photon-production events relative to the events re- 
sulting from ISR production, continuum production, and 
combinatoric background. 

Concerning the histogram of Fig. [3(c), the product 
r 77 (ry c (15)) x B{r) c (lS) -> K~ K + tt + tt~) is 1.7 ± 1.0 
times that for the Xco(l-P) state [lj], and in Fig. [3(c) the 
T] C (1S) signal contains ps 950 events (c/. the — 50 MeV/c 
interval of Fig. [3(a)), while the Xco(l-P) signal contains 
ss 550 events. It follows that the signal sizes agree well 
with the ratio expected on the basis of a two-photon 
production mechanism. In a similar vein, the ratio of 
the partial width T 11 (xcj) x B(y r ..i —> K~K + tt + it~) for 
Xco(l-P) and Xc2(lP) is 9±2 [lj], so that after taking into 
account the (2J+ 1) spin factors, the signals observed in 
Fig. [3(c) would be expected to be approximately in the 
ratio 1.8±0.4. The x C 2(l-P) signal contains w 200 events, 
and so is consistent with this expectation. 



IV. RECONSTRUCTION OF DD EVENTS 

Candidate DD events are reconstructed in the five 
combinations of D decay modes listed in Table Ql] (the 
use of charge conjugate states is implied throughout the 
text). Events are selected by requiring the exact number 
of charged-particle tracks defined by the relevant final 
state. 

Track selection requirements include transverse mo- 
mentum p t > 0.1 GeV/c, at least 12 coordinate measure- 



TABLE I: Charmonium states observed in the e + e — > 
K + K~n + n~X test data sample 





Mass [MeV/c 2 ] 






(2980.3 ±1.2) 


0-+ 


J/0 (IS) 


(3096.916 ± 0.011) 


1~ 


Xco(LP) 


(3414.75 ±0.31) 


0++ 


Xcl(lP) 


(3510.66 ±0.07) 


1++ 


Xc2(lP) 


(3556.20 ± 0.09) 


2++ 


Vc{2S) 


(3637 ± 4) 


o-+ 


4>(2S) 


(3686.09 ± 0.04) 


1"" 



TABLE II: D decay final states studied in this analysis; for 
channels N5, N6, and N7, inclusion of the corresponding 
charge conjugate combination is implied. 



Channel D decay mode D decay mode 



N4 


D"D" 


D u — > 


K--K+ 




■> A+7T- 


N5 


D°D° 


D° -» 


K~n+ 


D°- 


■> K+n-n° 


N6 


D°D° 


D° -» 


K--K+ 


D a - 


■> A + 7r _ 7r _ 7r + 


N7 


D°D° 


D° -> 


K~TT + TT + 7V 


D°- 


■> K+n-n a 


C6 


D+D- 


D+ - 


> K~1T + 1Y + 


D~ 


-)■ A + 7T _ 7r _ 


ments 


in the DCH, a 


maximum 


distance of closest ap- 



proach (DOCA) of 1.5 cm to the z-axis, with this point 
at a maximum DOCA of 10 cm to the xy-plane at z — 0. 

Kaon candidates are identified based on the normal- 
ized kaon, pion and proton likelihood values (Lk-, L v and 
L p ) obtained from the particle identification system, by 
requiring Lk/{Lk ± L-n) > 0.9 and Lk/{Lk ± L p ) > 
0.2. Tracks that fulfill L K /{L K + hj) < 0.82 and 
Lp/(L P ± L-x) < 0.98 are selected as pions. Addition- 
ally, in both cases the track should be inconsistent with 
electron identification. 

Photon candidates are selected when their deposited 
energy in the EMC is larger than 100 MeV. Neutral pions 
are reconstructed from pairs of photons with combined 
mass within [0.115,0.155] GeV/c 2 and a ir° mass con- 
straint is applied to them. 

The D candidate decay products are fitted to a com- 
mon vertex with a D meson mass constraint applied; can- 
didates with a x 2 fit probability greater than 0.1% are 
retained. Accepted DD pairs are refitted to a common 
vertex consistent with the e + e~ interaction region, and 
those with a \ 2 fit probability p v (D D) greater than 0.1% 
are retained. Events with ir° candidates other than those 
from a D or D decay of interest are rejected. These pre- 
selection criteria are identical for all five combinations of 
D decay modes. 

The signal regions for accepted, unconstrained D can- 
didates are then fitted using a multi-Gaussian signal func- 
tion 

rr*o 1 (m - mo) 2 
R{m) = / — -e ^~da (4) 
J<y ra 

with free parameters erg, r (minimal and maximal width) 
and mo; the background is described by a polynomial. 
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FIG. 2: (a) A'+A" 7r + 7r mass distribution for all events without any requirement on m 2 x ; (b) corresponding m x distribution; 
(c) m(K + K~n + n~) with the requirement m x < 10 (GeV/c 2 ) 2 ; (d) m(K + K~n + n~) with the requirement > 10 (GeV/c 2 ) 2 . 




Pi (K + KW)/GeV/c P( (K + KW) / GeV/c m^KW) [GeV/c 2 ] 



FIG. 3: Signal yield dependence on pt(K + K~ -k + tv~) for (a) ^(lS); (b) J/0. The vertical lines mark the pt(K + K~ tv + tt~) < 
0.05 MeV/c region for selecting two-photon events, (c) Resulting K + K~ tv + tt~ mass distribution after applying the principal 
selection criteria discussed in Sec. IIHI 



The full width at half maximum (FWHM) of the sig- 
nal lincshape in data is used to define each D signal re- 
gion; D candidates are selected from a region of width 
±1.5 FWHM around the mean mass. The mass windows 
are listed in Table Mil 

From the list of accepted DD candidates those pro- 
duced in two-photon events are then selected by apply- 
ing the three criteria defined in Sect. IIIII (summarized in 
Table HV]) . These criteria are also identical for all combi- 



nations of D decay modes. 

Depending on the decay mode, up to 2.5% of the events 
have multiple candidates which passed all selection crite- 
ria. In this case, the candidate with the best fit probabil- 
ity Pv(D D) is chosen. Based on MC studies, the correct 
candidate is selected in more than 99% of the cases with 
this method. 

The resulting invariant mass spectra for D meson can- 
didates after all selection criteria have been applied are 
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TABLE III: Summary of selection criteria used for identifying 
DD candidates. 



Channel 



number 

0f7T° 



D mass window D mass window 



[MeV/c 



[MeV/c 



N4 1863.4 ± 22 
N5 1 1863.4 ± 22 
N6 1863.4 ± 22 
N7 1 1863.4 ± 16 
C6 1868.5 ± 18 


1863.4 ± 22 
1863.4 ± 43 
1863.4 ± 16 

1863.4 ± 43 

1868.5 ± 18 


TABLE IV: Summary of requirements used for selecting only 
DD candidates from two-photon events. These criteria are 
identical for all decay modes. 


Channel wimiss 

[(GeV/c 2 ) 2 ] 


pt(DD) Eemc 
[GeV/c] [GeV] 


all modes > 10.0 


< 0.05 < 0.4 



as the experimental data. For each signal decay channel 
about 10 6 events were generated. Additional events were 
generated for background modes involving D* mesons. 
The GamGam two-photon event generator was used to sim- 
ulate 77 —> Z(3930) — > DD events, while the decays of 
the D and D mesons were generated by EvtGen [2(| ■ The 
detector response was simulated using the GEANT4 [2l| 
package. The program GamGam uses the BGMS formal- 
ism [23| . It was developed for CLEO and was used for 
example in the analysis of Xco(l-P): Xc2^P) — > 4-7T de- 
cays [23 ■ GamGam was later adapted to BABAR and used 
for the analysis of r) c (lS, 25) -> K^K^n* 

For small photon virtualities |qi,2| 2 (see Fig. [I} 
the differential cross section for the process e + e~ — > 
e + e~77, 77 — > X is given by the product Lx F x (7(77 — > 
X), where L is the two-photon flux. The form factor F 
extrapolates the process to virtual photons and is a priori 
not known. A plausible model 



shown in Fig. @] for events in which the mass of the re- 
coil D candidate lies within the defined signal region. In 
all modes, clear signals with small backgrounds are ob- 
tained. The resulting DD invariant mass distributions 
are shown in Fig. 03a) and EJb) for the neutral modes 
(N4, N5, N6, N7), and for the charged mode (C6), re- 
spectively. The combined spectrum is shown in Fig.[5{c). 
An enhancement near 3.93GeV/c 2 is visible. 

To estimate the amount of combinatoric background in 
the signal region, the two-dimensional space spanned by 
the invariant masses of the D and D candidates is divided 
into nine regions: one central signal region and eight side- 
band regions above and below the signal region as shown 
in Fig. El for the K~-k + -k + / K + -k~-k- (C6) mode. The 
mass range for the signal region is ± 1.5 FWHM around 
the mean mass. The sideband regions are 1.5 FWHM 
wide, leaving a gap of 1.5 FWHM between signal and 
sideband. No significant contribution from combinatoric 
background is observed in the DD spectrum (Fig. Efc)) . 

An attempt was made to isolate the signal in Fig. [5jc) 
by a weighting method. This assumes that signal and 
background events have different angular distributions, 
and was successfully used in a previous BABAR analy- 
sis [lj|. Simulations with a J PC = 2 ++ signal (gen- 
erated with its correct angular distribution) plus back- 
ground showed that the method works well with high 
signal statistics and moderate background, but is not re- 
liable with the limited statistics and background of the 
current analysis. Therefore, the method was not consid- 
ered further in the present analysis. 



V. MONTE CARLO STUDIES 



For modeling the detector resolution, efficiency stud- 
ies and the estimation of the two-photon width T 77 of 
the resonance, Monte Carlo (MC) events were generated 
which pass the same reconstruction and analysis chain 



F 



1 



1 - q{ /ml 



1 



1 - qi/ml 



(5) 



is used [2J|, with m v being the mass of an appropriate 
vector boson (p, J/tp, Z°). In the calculations relevant 
to this analysis m v = m(J/ip) was used, as the .£(3930) 
is expected to be a charmonium state. An alternative 
model was used in order to evaluate systematic uncer- 
tainties associated with MC simulations (see Section lXl) . 

To validate the GamGam generator its output was com- 
pared to that of another two-photon generator (TREPS) 
used by Belle 23 • The cross sections for the reac- 
■ e + e~77,77 — > r/ c (lS,2S) were calcu- 



tions e + e 



latcd in GamGam and compared to the Belle values [26| . 
In order to compare the different generators, the cross 
sections were calculated using the hypothetical values 
T 77 x B(r) c (lS,2S) -> final state) = 1 keV, and q\ 2 
was restricted to values smaller than 1 (GeV/c 2 ) 2 . The 
TREPS results were 2.11 pb for ^(15) and 0.86 pb for 
r] c (2S). The corresponding GamGam values were 2.13 pb 
and 0.84 pb, respectively. The two generators are in 
agreement at the level of a few percent. 

For a global check, the cross sections for the continuum 
reaction e + e~ — > e + e~77, 77 — > were calculated 

with GamGam for various CM energies and compared to 
QED predictio ns [ 27l . [28j , which describe the data with 
high accuracy [2{j]. Here, the agreement was slightly 
worse, due to the imperfect tuning of the GamGam pro- 
gram for these reactions. Similar results were obtained 
when checking against calculations with non-relativistic 
models for ?7 C (1S') and Xc2(l-P) [30| . Nevertheless this 
comparison showed that GamGam works properly under 
these conditions also. These studies lead to the assign- 
ment of a total systematic uncertainty of ±3 % associated 
with the MC simulation (see Sec. IXI[) . 
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FIG. 4: Candidate D invariant mass distributions after all selection criteria. The mass of the accompanying D candidate is 
required to lie within its signal region as defined in Table [TTTl (a) K~ir + in N4; (b) K + ir~ir° in N5; (c) K~ ir + ir + Tr~ in N6; (d) 
K+-K--K in N7; (e) K--k+ty + in C6. 
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FIG. 5: DD invariant mass distribution for the (a) D°D° and (b) D + D channels, (c) The combined DD invariant mass 
distribution. The open histogram in (c) shows the combinatoric background estimated from the D-mass sidebands. 



VI. PURITY OF THE 77 -> DD SAMPLE 

The selection criteria used to enhance the two-photon 
content of the DD sample were discussed in Sec. IIIII 
They were developed by investigating the reaction 
e+e" -> K+K-tt+tt'X of Eq. ©. Figure Ifc) shows 
that after the selection procedure the signals associated 
with 77-reactions, like that for the tj c (1S), are enhanced, 
while signals such as that for the J/ip , which are typical 



of ISR production, are suppressed. The pt(DD) distri- 
bution is shown in Fig.[7]for events in the Z(3930) signal 
region, defined as the region from 3.91 to 3.95GeV/c 2 . 
Here the p t (DD) selection criterion has not been applied. 
The data are fitted with a curve for 77 events obtained 
from MC, plus a linear background derived from sideband 
studies of the DD mass spectrum. The fit indicates that 
the majority of DD candidates in the signal region result 
from two-photon interactions. 
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of Table [Til The boxes correspond to the signal and sideband 
regions. 



VII. RECONSTRUCTION EFFICIENCY 

The reconstruction efficiency for each decay mode is 
calculated function of m(DD) using MC events 

which pass the same reconstruction and selection cri- 
teria as real events and includes detector acceptance, 
track reconstruction- and particle identification efficien- 
cies. The mass-dependent efficiency €i(m(DD)) for each 
channel i is fitted with a polynomial in m(DD) and 
is found in each case to decrease with increasing DD 
mass. For the combination of modes (Fig.[5jc)), an over- 
all weighted efficiency e B (m(DD)), which includes the 
branching fractions for the D decays, is computed using 



\m(DD)) 



5 U=iNi(rn(DD)) 
2 



E 



5 Nj(m(DD)) 
1 ef (m(DD)) 



(0) 




0.1 0.15 0.2 
p (DD) / GeV/c 



FIG. 7: Distribution of p t (DD) for data in the Z(3930) sig- 
nal region. The fitted lineshape consists of the expected 77 
lineshape obtained from MC plus a linear background (dotted 
line) . The vertical line shows the pt criterion for selecting 77 
events. The histogram shows the shape of the pt(DD) distri- 
bution from simulated D* D events with missing 7r° or 7. The 
bump in this distribution is not seen in the data distribution, 
indicating that any D* D background is small. 
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FIG. 8: Mass dependence of the weighted reconstruction effi- 
ciency calculated using Eq. ©, described by a straight line. 



as was done in Ref. [3l|; Ni(m(DD)) is the number of 
DD candidates in the data mass spectrum for channel 
i, and ef (m(DD)) is denned as the product of the effi- 
ciency €i as parameterized by the fitted polynomial and 
the branching fraction Bi for the i-th channel, as follows 



\m{DD)) = ei(m(DD)) x B % 



(7) 



The factor \ originates from referring to DD (D°D° 
and D + D~) events; the factor 5 from summing over the 
five channels. Figure |S] shows the mass dependence of 
e B (m(DD)), which is parameterized by a straight line. 
The large uncertainties are due to the limited statistics 
available in the data samples. The error bars do not con- 
tain the uncertainties in the branching fractions; these 
will be discussed separately in Sec. I XII in the context of 
systematic error estimation. The data are weighted by 
this mean efficiency, which is scaled by a constant value 
d to obtain weights near one, 



e(m(DD)) =dx e B (m{DD)) 



(8) 



as weights far from one might result in incorrect errors 
for the signal yield obtained in the maximum likelihood 
fit [H]. The resulting DD mass distribution will be dis- 
cussed in Sec. IVIIII 



VIII. DETECTOR RESOLUTION AND SIGNAL 
YIELD 

Monte Carlo events are used for the calculation of the 
mass-dependent detector resolution. The mass resolu- 
tion is determined by studying the difference between 
the reconstructed and the generated DD mass (Am res ). 
As an example, the distribution for channel C6 is shown 
in Fig. inja). A good description of the distribution is 
obtained using a multi-Gaussian fit (Eq. (|4])). The pa- 
rameters r and a (m(DD)) were determined for every 
decay channel. The variation of <ro(m(DD)), which is 
parameterized by a second order polynomial, and of the 
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width (FWHM) of the resolution function with increasing 
mass are shown in Figs. HJb) and^c). For channel C6, 
r = 5.380±0.137 and a (m(DD))_ = (-0.038 + 0.018to- 
0.002m 2 ) GcV/c 2 , where m(DD) is given in units of 
GeV/c 2 . The distributions of Fig. [9] are well-described 
by the fitted curves shown. Comparing the generated 
Z(3930) mass with the reconstructed MC value shows 
that the latter is systematically low by about 0.9 MeV/c 2 , 
independently of the fit model. This effect is observed 
both in the combined fit and in fits to the individual 
channels. The measured J/ip mass in the K + K~tt + it~ 
test sample (Sect. [nB differs by the same value from 
the world average [PM; t his offset has been seen in other 
77 studies at BABAR jig ] as well. Accordingly, the mass 
value obtained from the fit to data will be corrected by 
+0.9 McV/c 2 . This offset value will also be used as a con- 
servative estimate of the systematic uncertainty in the 
mass scale. The difference between the generated and 
reconstructed decay width values amounts to 0.14 MeV, 
and is discussed in Sec. I XII with respect to systematic 
error estimation. 

In order to describe the signal structure in data 
around 3.93GcV/c 2 a relativistic Breit-Wigner function 
BW(m) is used, where 



BWirn) 





, 2L+1 


( Pm S 


) © 


\Pmo , 





F 



V to / (to 2 , — to 2 ) 2 + r 2 



rn'"0 



(9) 

with too as the nominal mass of the resonance; the Blatt- 
Weisskopf coefficients F r for different angular momentum 
values L are given by 



F r (L = 0) 
F r (L = 1) 

F r (L = 2) 

and the value 



1 



y/l + {Rpra? 

v/9 + 3(i?p mo ) 



(RPr, 



^9 + 3(R Pm ) 2 + (Rp„ 



(10) 

(11) 

(12) 



as discussed previously in this section. The background 
is parameterized by the function 

D(m) oc \J to 2 — to 2 (m — m t ) a exp [— /3(m — m t )] (14) 

which takes the DD threshold m t into account. In the 
lower mass region, the lineshape does not describe the 
background exactly. Other functional forms were tried 
fSec. IXl) ). but no improvement was obtained. The data 
and the curves which result from the standard (J = 2) 
fit are shown in Fig. [TU1 

From the unbinned maximum likelihood fit to the 
five mass spectra the Z(3930) values too = (3925.8 ± 
2.7) MeV/c 2 and T r = (21.3 ± 6.8) MeV are obtained for 
the mass and total width, respectively (all errors in this 
section are statistical only). The mass is corrected by 
+0.9 MeV/ c 2 as described above, resulting in a final mass 
value of (3926.7 ±2.7) MeV/c 2 . The efficiency-corrected 
yield amounts to TV = (76+ 17) signal events. This value 
is based on weights around 1 as discussed in Sec. IVIII 
taking the constant used to scale the efficiency into ac- 
count (see Eq. [5]), this corresponds to a total Z(3930) 
signal of N e s = (285 ± 64) x 10 3 events. 

The statistical significance of the peak is 5.8a and is 
derived from the difference A In £ between the negative 
logarithmic likelihood of the nominal fit and that of a fit 
where the parameter for the signal yield is fixed to zero. 
This is then used to evaluate a p-value: 



P = I f(z;n d )dz 

/2Aln£ 



(15) 



where f(z; rid) is the x 2 PDF and rid is the number of 
degrees of freedom, three in this case. We then determine 
the equivalent one-dimensional significance from this p- 
value. 



IX. ANGULAR DISTRIBUTION AND SPIN OF 
THE Z(3930) STATE 



R=1.5 (GeV/c)" 1 

is used, corresponding to the value given in Ref. (33j . The 
mass-dependent width is given by 



/ \ 2L+1 , 

( HHL) ( m ° 

\Pmo 



F 



(13) 



with r. r the total width of the resonance. Here the ex- 
istence of other possible decay modes is ignored. The 
momentum of a given D candidate in the DD center of 
mass frame is denoted by p m ; p mo is the corresponding 
value for to = too. In the standard fit, spin J = 2 (L = 2) 
is chosen on the basis of the angular distribution analysis 
described in Sec. IIXI 

The signal function is convolved with the mass- and 
decay-mode-dependent resolution model parameterized 



General conservation laws limit the possibilities for the 
J PC values of the Z(3930) state. For two-photon produc- 
tion the initial state has positive C-parity and hence the 
final state must have positive C-parity also. For the DD 
final state, C = (— 1) L+S = (-1) L since the total spin 
S is zero. Positive C-parity then implies that the DD 
system must have orbital angular momentum L which is 
even, and hence have even parity. It follows that for the 
Z(3930) state J PC = J++ with J = 0,2,4... In order 
to investigate the possible values of J, we have compared 
the decay angular distribution measured in the Z(3930) 
signal region to the distributions expected for J = and 
J = 2; higher spin values are very unlikely for a state 
only 200 MeV/c 2 above threshold. 

The decay angle 9 is defined as the angle of the D 
meson in the DD system relative to the DD lab momen- 
tum vector. Figure [TT] shows the if (3930) signal yield 
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FIG. 9: (a) Detector resolution Am res for channel C6; the fitted curve is described in the text; (b) mass-dependence of the 
resolution function parameter <tq; (c) mass-dependence of the FWHM of the resolution function. 




4 4.2 
m(DD) [GeV/c 2 ] 



FIG. 10: Efficiency-corrected mean DD mass distribution 
with standard fit. The dashed curve shows the background 
lineshape (see Sec. IVIH"]) . 



obtained from fits to the DD mass spectrum for ten re- 
gions of |cos0|. The data have been weighted by a cos0- 
dependent efficiency, which was determined in a similar 
manner as described in Sec. IVIll for the mass-dependent 
efficiency (Fig. [T2]) . In these fits, the mass and width 
of the resonance have been fixed to the values found in 
Sec. I Villi and Eq. (|T4f has been used to describe the 
background. Other background models have been tried, 
obtaining distributions fully consistent with Fig. [IT] 

The function describing the decay angular distribution 
for spin 2 has been calculated using the hclicity formalism 
and has the form 




0.2 0.4 0.6 0.8 1 
IcosOl 



FIG. 11: Signal yield as a function of |cos#| derived from fits 
to the efficiency-corrected DD spectrum. The solid curve is 
the expected distribution for spin 2 with dominating helicity-2 
contribution, the dotted straight line is for spin 0. 



It has been assumed that the dominating amplitude has 
hclicity 2. This is in agreement with previous measure- 
ments [HI and theoretical predictions [24| H(| ■ The dis- 
tribution of Eq. (|16|) was fitted to the experimental an- 
gular distribution, and a x 2 /NDF value of 5.63/9 was 
obtained, with NDF indicating the number of degrees of 
freedom. For a flat distribution, which is expected for 
spin 0, a x 2 /NDF = 15.55/9 was obtained. It follows 
that the preferred J PC assignment is 2 ++ . 



X. TWO-PHOTON WIDTH OF THE Z(3930) 
STATE 



From the efficiency-corrected number of observed sig- 
nal events, N e B , we determine the total experimental 
cross-section 



cr e xp(e + e" 



dN 
dcos I 



cx sin 8. 



(16) 



■e+e~77,77 -> Z(3930), Z(3930) -> DD) 
N eB / J Cdt = 741±166fb (17) 
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FIG. 12: Angular dependence of the weighted reconstruction 
efficiency e s (cos#) based on Eq. \E[ described by a second 
order polynomial. 



where the integrated luminosity for the data sample an- 
alyzed is J Cdt = (384 ± 4) fb _1 and the error is only 
statistical. 

On the other hand, the cross-section for Z(3930) pro- 
duction is given by 



with 



77,77^ Z(3930)) = LxFxa(jj^ Z(3930)) 

(18) 



cr(77 -> Z(3930)) = / 4tt(2J + l)(7ic) 2 10- 6 m| x 



r-. 



r, 



Km (m 2 z 



m 2 ) 2 +m 2 T 2 ot 



dm 2 (19) 



and can be calculated using GamGam. Here L is the two- 
photon flux, F is the form factor (see Sec.|V|, mz (r to t) 
is the resonance mass (width), and T 77 is the two-photon 
width of the resonance. The kincmatical factor K is given 
by K = {q\q2) 2 — q 2 q 2 {qi represent four vectors of pho- 
tons). Further information can be found in Rcfs. [351 ] . 
[22j and [24]]. The cross-section depends on the spin of 
the resonance and on T 77 . It is plotted for J = 2 and 
J = in Fig. [13] as a function of r 77 . From a com- 
parison to the experimental cross-section (Eq. (fTT)) '). the 
partial width T 77 x S(Z(3930) -> DD) is found to have 
the value (0.24 ± 0.05) keV when J = 2 is chosen as the 
most probable spin value (see Sec. IIX]) . 



XI. SYSTEMATIC ERROR ESTIMATION 

Several sources of systematic uncertainty have been 
considered for the mass, decay width, and signal yield 
of the Z(3930) state. The yield determines the value 
of T 77 x £(Z(3930) -> DD). The standard fit to the 
efficiency-corrected mass spectrum is repeated with ap- 
propriate modifications. The differences A between the 
results obtained and the standard results arc used as es- 
timates of systematic uncertainty. No correlations have 
been taken into account. The results are summarized in 



FIG. 13: Dependence of the cross-section a(e + e — » 77 — > 
Z(3930) -> DD) on the two-photon width T 77 x £(Z(3930) ->• 
DD), calculated with the two-photon generator GamGam. The 
upper solid line is for spin 2, while the lower solid line is for 
spin 0. The measured value (horizontal dashed line) and its 
uncertainty range (horizontal dot-dashed lines) are indicated. 



Table [V] Deviations for the mass (|Am|), total width 
(|Ar|) and two-photon width (|A(r 77 x B)\) are consid- 
ered negligible if they are less than 0.05 MeV/c 2 , 0.05 McV 
and 0.0005 keV, respectively. 



A. Fit parameterization 

Signal Lineshape: The standard fit has assumed 
spin J = 2 for the resonance (Sec. IVIII]) . Using different 
spin values and R values has no significant impact on the 
results (Table [V] A(r 77 x B) numbers are given for spin 
,7 = 2 only). 

Background Description: Different paramcteriza- 
tions of the background in the m(DD) distribution have 
been used. Besides the nominal background (Eq. (fT4]l). 
the following background shape was tried 



D' '(m) oc I 1 — exp 







1 m 
7 - -1 
\ a 



(20) 

the fit had a slightly worse, but still acceptable, likelihood 
value. The mass value changes by Am = +0.4MeV/c 2 , 
the width by Ar = +3.0 MeV, the signal yield by +9 
events with respect to the standard fit, and T 77 x B 
chang es accordingly by +0.029 keV (Table |V|. Other 
background models yield consistent estimates for this 
source of systematic uncertainty. 



B. Detector resolution 

Fit Precision and Mass Scale: A fit of the convo- 
lution of signal lineshape and resolution model to the MC 
sample has been performed. The mass offset observed in 
MC has been included by correcting the mass value by 
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+0.9 McV/c 2 . As a conservative estimate, this number 
is also used as the systematic uncertainty for the mass 
scale. The deviation between the generated width and 
the value obtained from the fit is 0.14 MeV, and again 
this is used as a conservative estimate of systematic un- 
certainty. Based on the uncertainty of the width, a value 
of Ar 77 xg = 0.001 keV is derived. 

Resolution Model: The parameters of the multi- 
Gaussian resolution model were modified. The number 
of steps was enlarged from 25 to 35, the total convolu- 
tion range for each data point enlarged by +0.02 MeV/c 2 , 
and the parameter r of the multi-Gaussian was varied 
within its fit uncertainty Sr. The corresponding shifts in 
the mass are Am = +0.2, < 0.05 and < 0.05McV/c 2 . 
For AT, shifts of -0.2, -0.9 and -0.1 MeV are obtained; 
from the modified signal yield, shifts of —0.003, —0.003 
and < 0.0005 keV were obtained for T 77 x B (Table EJ- 

C. Combined reconstruction efficiency 

Parameterization: The average mass-dependent 
reconstruction efficiency has been parameterized by a 
straight line in the standard fit (Fig. [5]). Using a fit 
with a second order polynomial, the width changes by 
—0.4 MeV; no mass shift was observed with respect to 
the standard fit result. For the signal yield, +1 entry 
is obtained; this yields no significant shift for r 77 x B 
(Table E| . 

Tracking and Neutrals Correction: For the 

tracking efficiency a correction by —0.8% is applied per 
charged-particle track. This gives a correction factor of 
0.968 for modes N4, N5, and 0.953 for N6, N7 and C6. 
The systematic uncertainty assigned to the tracking ef- 
ficiency is 1.4% per track for decays with more than 5 
charged particle tracks and 1.3% otherwise. The result- 
ing uncertainty for T 77 x B is 0.022 keV. Concerning effi- 
ciency corrections for neutral particles, a correction fac- 
tor of 0.984 with an uncertainty of 3% per 7r° is used for 
modes N5 and N7. The resulting uncertainty for T 77 x B 
is 0.003 keV (Table El- 
Uncertainty on the D Branching Fractions: 
The errors on the D branching fractions have been taken 
into accout by varying the values of Bi used in Eq. ^ 
within their standard deviations. No significant change 
is observed in mass and decay width. For the two-photon 
width A(r 77 x B) = ±0.010 keV is obtained (Table El- 
Effect of Angular Distribution on Efficiency: 
The MC data sample used to obtain the efficiency and 
resolution was generated with a flat distribution in cos 9. 
To estimate the effect of the angular distribution on the 
reconstruction efficiency, a MC sample described by a 
sin 4 9 distribution has been generated and reconstructed. 



Comparing these reconstructed data with the nominal 
MC sample, the mean efficiencies differ by 8%, relatively, 
resulting in A(r 77 xB) = ±0.018 keV. 

D. Cross-section calculation from GamGam 

Precision: In Sec. El a relative uncertainty of ±3 % 
was obtained for the calculated cross-section. Propagat- 
ing this error into the calculation of T 77 x B, an uncer- 
tainty A(r 77 X B) = ±0.007 keV results. 

Form Factor: In the standard analysis the form fac- 
tor of Eq. ([5]) has been used with m v = m(J/ip). In order 
to estimate potential systematic effects, the cross-section 
was evaluated using a model predicted by perturbativc 

QCD m 



(1 - ql/ml - ql/mlf 

The cross-section calculated with GamGam does not in- 
crease significantly (sa 0.1%) compared to that obtained 
using Eq. ([5]). Simultaneously the experimental efficiency 
decreases by 1%, so that the net effect on T 77 x B is small. 
Similar effects have been observed when data and calcu- 
lations with and without q 2 selection criteria are com- 
pared [25|, [26|, [3?J , and also in a previous CLEO analy- 
sis (38[. As a result a systematic uncertainty of ±1% is 
attributed to form factor uncertainty and this yields a 
deviation A(r 77 x B) = ±0.002 keV 

E. Other uncertainties 

Particle Identification (PID): For PID studies, 
the pion selection criteria have been tightened signifi- 
cantly, and the efficiency has been recalculated accord- 
ingly. The fit to the mass spectrum yields a change of 
-0.4MeV/c 2 for the mass and -1.8 MeV for the width. 
For A(r 77 x B) a, change of -0.004 keV results. 

D Mass Uncertainty: The uncertainty of the D 
meson mass is taken into account. Both for D° and D^, 
the uncertainty is 0.17MeV/c 2 [l4[, which results in an 
uncertainty of ±0.34MeV/c 2 in the mass of the Z(3930) 
state. 

Integrated Luminosity Uncertainty: For the in- 
tegrated luminosity, an uncertainty of ±1% is assigned. 
From this an uncertainty A(r 77 X B) = ±0.002 keV is 
obtained. 



F. Total systematic uncertainty 

The systematic uncertainty estimates discussed in 
Sees. IXI A1 - IXI El arc summarized in Table El The indi- 
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TABLE V: Results of the systematic uncertainty studies for the mass, decay width and efficiency-corrected signal yield of the 
Z(3930) state. Listed are the differences with respect to the standard values. A(F 77 x B) numbers are given for spin J = 2 
only. For the combined error, the values are added in quadrature. 



Source of Systematic Uncertainty 


Ara(Z(3930)) 


Ar(Z(3930)) 


A(r 77 x B) 




[MeV/c 2 ] 


[MeV] 


[keV] 


Choice of Spin J = 1, J = 


< 0.05 


< 0.05 


- 


Value of R (Breit-Wigner) 
Background D'(m) 


< 0.05 


< 0.05 


< 0.0005 


0.4 


3.0 


0.029 


Fit precision and mass scale 


0.9 


0.1 


0.001 


Convolution steps = 35 


0.2 


0.2 


0.003 


Convolution range +0.02 MeV/c 2 


< 0.05 


0.9 


0.003 


Resolution multi- Gauss r ± Sr 


< 0.05 


0.1 


< 0.0005 


Combined reconstr. efficiency: polynomial 


< 0.05 


0.4 


< 0.0005 


Tracking efficiency correction 


< 0.05 


< 0.05 


0.022 


7r^ efficiency correction 


< 0.05 


< 0.05 


0.003 


Error in D branching fractions 


< 0.05 


< 0.05 


0.010 


Efficiency: angular distribution 






0.018 


Generator precision 






0.007 


Choice of Form Factor 






0.002 


P1D 


0.4 


1.8 


0.004 


Uncertainty in D mass 


0.3 






Luminosity 






0.002 


Combined error 


±1.1 


±3.6 


±0.04 



vidual estimates are combined in quadrature to yield net 
systematic uncertainty estimates on the Z(3930) mass, 
total width and value of T 77 x S(Z(3930) DD) of 
l.lMeV/c 2 , 3.6 MeV and 0.04 keV, respectively, as re- 
ported on the last line of Table [V] 

XII. SUMMARY 

In the 77 — > DD reaction a signal in the DD mass 
spectrum has been observed near 3.93 GeV/c 2 with a sig- 
nificance of 5.8cr which agrees with the observation of 
the Z(3930) resonance by the Belle Collaboration (l3| . 
The mass and total width of the Z(3930) state are mea- 
sured to be (3926.7 ± 2. 7(stat) ± l.l(syst)) MeV/c 2 and 
(21.3 ± 6.8(stat) ± 3.6(syst)) MeV, respectively. 

The production and decay mechanisms allow only pos- 
itive parity and C-parity, and an analysis of the Z(3930) 
decay angular distribution favors a tensor over a scalar 
interpretation. The preferred assignment for spin and 
parity of the Z(3930) state is therefore J p _^ = 2++. 
The product of the branching fraction to DD times the 
two-photon width of the Z(3930) state is measured to 
be T 77 x £(Z(3930) ->• DD) = (0.24 ± 0.05(stat) ± 
0.04(syst)) keV, assuming spin J = 2. The parameters 
obtained are consistenwith the Belle results, and with 
the expectations for the Xc2(2P) state. 
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